clc;
clear;

syms x y
u = -(2-pi*sin(pi*x)).*cos(2*pi*y);
u_x = diff(u, x);
u_y = diff(u, y);
u_x = simplify(u_x);
u_y = simplify(u_y);
disp(u_x);
disp(u_y);